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1 2-Chemokine Gene Signature Identifies 
Lymph Node-like Structures in 
Melanoma: Potential for Patient Selection 
for Immunotherapy? 

Jane L. Messina' *, David A. Fenstermacher^, Steven Eschrich^, XiaotaoQu^, Anders E. Berglund^, 
Mark C. Lloyd^ Michael J. SchelP, Vernon K. Sondak' Jeffrey 5. Weber' & James J. Mule' -^ 

^ Cutaneous Oncology, "^Biomedical Informatics, ^Biostatistics, "^Immunology, ^the Analytic Microscopy Core, Moffitt Cancer Center, 
Tampa, PL, USA, ^the Departments of Pathology and Cell Biology, University of South Florida, Tampa, FL, USA. 

We have interrogated a 12-chemokine gene expression signature (GES) on genomic arrays of 14,492 distinct 
solid tumors and show broad distribution across different histologies. We hypothesized that this 
12-chemokine GES might accurately predict a unique intratumoral immune reaction in stage IV 
(non-locoregional) melanoma metastases. The 12-chemokine GES predicted the presence of unique, lymph 
node-like structures, containing CD20"^ B cell follicles with prominent areas of CD 3^ T cells (both CD4^ and 
CD8^ subsets). CD86^, but not FoxP3^, cells were present within these unique structures as well. The direct 
correlation between the 12-chemoldne GES score and the presence of unique, lymph nodal structures was 
also associated with better overall survival of the subset of melanoma patients. The use of this novel 
12-chemoldne GES may reveal basic information on in situ mechanisms of the anti- tumor immune 
response, potentially leading to improvements in the identification and selection of melanoma patients most 
suitable for immunotherapy. 

Skin cancers are the most common of all human malignancies. Although melanoma constitutes a minority of 
skin cancers, it is considered to be the most deadly. The treatment of primary melanoma is predominantly 
surgical. Once melanoma has metastasized to lymph nodes, radical surgery to remove the nodes in the 
affected basin is performed, often followed by immunotherapy and sometimes by radiation. Melanoma metastatic 
to organ sites beyond the regional lymph nodes requires, in most cases, systemic therapy. Melanoma is notor- 
iously insensitive to standard chemotherapy drugs that are widely used for other forms of cancer. However, 
melanoma is occasionally, and sometimes dramatically, responsive to immunotherapy. 

Immunotherapy is fast becoming an important part of the treatment armamentarium for advanced melanoma, 
but the degree of its clinical effectiveness varies among patients. In spite of the clinical success observed with 
antibodies against CTLA-4 and cytokines (e.g., high-dose IL-2^), as well as the adoptive transfer of tumor 

infiltrating lymphocytes^'^, dramatic responses are observed in a minority of patients, while the majority of 
patients treated with those agents do not show significant clinical benefit. This current limitation underscores 
the need for the discovery of immune-related biomarkers and gene expression signatures that can identify/predict 
the likelihood of melanoma patients achieving a robust response of prolonged duration to immunotherapy. 

It is becoming widely accepted that immune cell infiltrates in human melanoma and other solid tumors have 
prognostic value^'^. Many of these studies described the subset composition of the lymphocytic infiltrates that 
were either distributed diffusely within tumor parenchyma or localized diffusely to the peripheral rim of the mass 
and the tumor/stroma interface. In addition, recent attention has been focused on defining and potentially 
utilizing a so-called "immune score" based on the types and degrees of immune cell infiltrates for classification 
of human cancer^. 

We have been studying the intriguing observation that ectopic lymph node-like structures reside within the 
parenchyma of some solid tumor masses^'^°. Our initial work focused on invasive colorectal carcinoma (CRC)^\ 
We found that a metagene (i.e. several tightly correlated sets of genes) score for inflammation and immunity 
accurately predicted host immune reaction when applied to this tumor type. Indeed, we identified, through the 
use of a 12-chemokine GES, a subset of CRCs that contained unique, tumor-localized ectopic lymph node-like 
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structures (TL-ELNs). This molecular classifier and the presence of 
these lymphoid structures were associated with better patient out- 
come, irrespective of other defined clinical and biologic variables. 
Based on their known selective capacity to potently recruit and con- 
centrate immune cell subsets in vivo, and also their known asso- 
ciation with peripheral lymph node formation and function^^"^^, it 
is possible that chemokines locally produced within the micro- 
environment of the tumor mass are responsible for the formation 
of these ectopic lymph node -like structures. It is also conceivable that 
these lymph nodal structures are active and are playing an important, 
positive role in eliciting an endogenous anti-tumor immune response 
(initially local but becoming systemic) that results in improved 
patient survival. 

In the current study, the 12-chemokine GES has been interrogated 
on genomic arrays of 14,492 distinct solid tumors (primaries and 
metastases) and is shown to have broad distribution across different 
histologies. Further, we show that this signature can, in fact, identify 
the presence of unique, TL-ELNs in melanoma metastases, which 
also appears to associate with better patient outcome. The implica- 
tion of these findings for current and future immunotherapies for 
melanoma is discussed. 

Results 

We had earlier identified in colorectal cancer a unique 12-chemokine 
(CCL2, CCL3, CCL4, CCL5, CCL8, CCL18, CCL19, CCL21, CXCL9, 
CXCLIO, CXCLll, and CXCL13) GES from a metagene grouping 
with overwhelming enrichment for immune-related and inflam- 
mation-related genes^\ In the current work, we interrogated this 
signature on 14,492 distinct solid tumors (primaries and metastases) 
with at least 30 per tumor type using gene chip technology. For this 
study, we employed principal component analysis (PCA), calculating 
a PCA score based on the first principal component - PCI - that 
approximates an average Z-score based on the log2 intensities for 
each of the 12 genes in the chemokine GES. Figure 1 shows the 
distribution of the 12-chemokine GES across different solid tumor 
histologies. 

Gene expression for the 12 chemokine genes was combined using 
PCA to generate a chemokine signature score for 120 individual, 
stage IV patient non-locoregional metastasis samples (Figure 2a.). 
The first principal component (—57% variance) is used to represent 
the chemokine signature. Melanomas with a high chemokine sig- 
nature showed increased expression values in all 12 chemokine genes 
when compared to melanomas with a low chemokine signature. 
Statistical tests of differences among the 12 genes confirmed the 
PCA results, indicating up-regulation of gene expression (p value 
ranges, 3.25e-05 to 2.22e-16, Figure 2b.). Melanomas that had the 
highest (n=9; red) or lowest (n=12; green) chemokine signatures 
were identified (Figure 2c.). 

Microscopically, representative tissue sections of resected meta- 
static melanoma specimens from patients with the lowest 12-chemo- 
kine gene scores revealed, in all cases, a minimally to absent 
lymphocytic peritumoral host response, with low to no appreciable 
expression of lymphocytic markers (not shown). Conversely, the 
highest 12-chemokine gene-scoring melanoma samples revealed, 
in all cases, a marked peritumoral lymphocytic host response con- 
sistent with ectopic lymph node-like structures, evident by H&E 
staining (Figure 3a.,b.,c.) and by IHC (Figure 4 a.,b.,c.,d.), particu- 
larly at the interface between tumor and stroma. Ectopic lymph 
node-like structures were observed intratumorally as well, occasion- 
ally accompanied by a low number of tumor infiltrating lympho- 
cytes. The tumor-localized ectopic lymph node-like structures (TL- 
ELNs) were found to contain prominent lymphoid follicles. As 
shown in Figure 4, CD20^ B cells were present almost exclusively 
within the follicular structures. The majority of the CD 3^ T cells 
(both CD4^ and CD8^ subsets) were concentrated in parafolUcular 
cortex-like zones; scattered T cells were also seen within these 



Tissue Type 


Total # 
CEL files 


# Above 
90th 

percentile 


0/ AK^oA 

/o Above 
90th 

percentile 


exact_test p 
value 


Oral Cavity 


98 


25 


25.51 


8.08E-06 


Cervix 


75 


19 


25.33 


0.000106163 


Tongue 


32 


8 


25 


0.011540008 


Skin* 


569 


115 


20.21 


1.11E-13 


Lung 


2708 


488 


18.02 


3.63E-47 


Soft Tissue 


97 


14 


14.43 


0.170421294 


Bladder 


212 


27 


12.74 


0.20279835 


Larynx 


56 


7 


12.5 


0.501205795 


Breast 


3705 


401 


10.82 


0.052757354 


Stomach 


133 


13 


9.77 


1 


Large Bowel 


2111 


169 


8.01 


0.000837406 


Kidney 


850 


61 


7.18 


0.003839317 


Thyroid 


71 


5 


7.04 


0.550835068 


Esophagus 


90 


6 


6.67 


0.377460472 


Rectum -An us 


188 


10 


5.32 


0.027468551 


Liver 


116 


5 


4.31 


0.041696746 


Endometrium 


333 


14 


4.2 


0.000131198 


Small Intestine 


52 


2 


3.85 


0.167045365 


Uterus 


377 


13 


3.45 


2.00E-06 


Pancreas 


468 


15 


3.21 


2.70E-08 


Ovary 


670 


21 


3.13 


8.78E-12 


Renal Pelvis 


62 


1 


1.61 


0.01881245 


Brain 


438 


4 


0.91 


2.30E-15 


Prostate** 


981 


5 


0.51 


9.64E-39 


Total 


14492 


1448 







Figure 1 | Interrogation of a 12-chemokine gene expression signature 
across solid tumors of differing histology: Principal component analysis 
(PCA). Gene expression microarrays on 14,492 samples (primaries and 
metastases) were examined and total PCA score was derived as a measure 
of chemokine signal. The top 10% of samples were selected as high. For this 
study, we calculated a PCA score based on PCI that approximates an 
average Z-score based on the log2 intensities for each of the 12 chemokine 
genes in the gene expression signature. RNA was extracted from tumors 
frozen within 15 minutes of surgical removal (according to strict standard 
operating procedures). *The 569 skin tumor samples include non- 
melanomas (i.e. basal cell, squamous cell, and Merkel cell carcinomas), 
primary melanomas, as well as 120 non-locoregional metastases from stage 
IV patients. Compared to all other samples, RNA was extracted from 
paraffin of suboptimal amounts of tumor tissue. 

tumors. CD86^ antigen presenting/dendritic cells were present 
within the structures as well; however, FoxP3^ cells (putative T reg- 
ulatory cells) were found to be very few in numbers. 

Clinical information was available on 10 of the 21 selected mela- 
noma patients. Table 1 shows the overall survival of these patients 
and the 12-chemokine GES that identifies the presence or absence of 
TL-ELNs for these 10 cases. There was a highly significant (p = 0.008) 
association between increased survival and the value of the mean 
score of the 12-chemokine GES. We also obtained descriptive 
information between the score of the 12-chemokine GES, site of 
the metastasis, and treatment received, which did not appear to show 
associations. The patients were those who predominantly had their 
metastatic lesions resected to render them free of disease, and for 
whom long survival can often be seen, especially from the time of the 
primary diagnosis and death. Of interest, one patient with a long 
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Principal Component Analysis of Melanoma Patients 
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Gene 


Welch T-test 
P-value 


GSE LOW 


GSE High 


Mean 


STD 


Mean 


STD 


CCL18 


3.25E-05 


6.02 


2.33 


10.78 


1.42 


CCL19 


4.42E-11 


4.85 


0.48 


9.49 


1.08 


CCL2 


1.58E-06 


7.92 


1.27 


11.41 


0.98 


CCL21 


1.38E-05 


6.05 


1.06 


10.14 


2.13 


CCL3 


5.12E-07 


5.18 


1.20 


8.80 


0.98 


CCL4 


1.03E-09 


7.22 


1.09 


11.62 


0.55 


CCL5 


4.36E-13 


6.90 


0.68 


11.72 


0.56 


CCL8 


l.lOE-11 


6.14 


0.92 


11.68 


0.80 


CXCLIO 


1.46E-11 


6.15 


1.23 


12.24 


0.43 


CXCLll 


2.22E-16 


3.31 


0.25 


8.73 


0.66 


CXCL11_2 


2.09E-14 


4.42 


0.62 


10.21 


0.67 


CXCL13 


4.73E-14 


4.40 


0.61 


10.71 


0.87 


CXCL9 


4.10E-10 


5.71 


1.55 


12.25 


0.73 



c. 




Figure 2 | Melanoma metastasis gene expression of 12 chemokines. 

(a.) Gene expression for the 12 chemokine genes ( 13 Affymetrix probe sets) 
was combined using PGA. The scatter plot consists of melanoma samples 
plotted by the first two principal components. The first component, 
representing 57% of the variability within samples, was used as the 
chemokine signature score (GSE) for 120 non-locoregional melanoma 
metastases samples, (b.) The melanoma samples with high chemokine 
signature scores (n=9) showed significantly increased expression values in 
all 12 chemokine genes (13 probe sets) when compared to melanomas with 
low chemokine signature (n= 12), as shown by mean and standard 
deviation of each gene, together with p value for differences, (c.) Patients 
whose melanoma samples had the very highest (n=9, red) or lowest 
(n= 12, green) chemokine signatures are shown in the heatmap with the 
expression of the 12 chemokine genes (13 probe sets). Samples are labeled 
by rank order; samples 1-12 have low chemokine GSE, and samples 112- 
120 have high GSE. 




Figure 3 | Presence of ectopic lymph node-like structures in melanoma. 

The nine highest 12 -chemokine gene expression signature scored 
melanomas revealed a marked peritumoral lymphocytic host response, 
organized as ectopic lymph node-like structures by H&E staining (a.,b.,c.). 
Arrows (a. and b.) depict individual ectopic lymph node-like structures; 
Gircle (b.) depicts an area of numerous ectopic lymph node-like structures. 

survival duration had received immunotherapy with an anti-CTLA4 
antibody (ipilimumab) and has had a partial response of over 30 
months duration. 

Discussion 

We have hypothesized that immune gene-related expression signa- 
tures can predict the presence of unique histological features of 
lymphoid cell infiltrates in solid tumor masses that correlate with 
better overall patient survival. In the current pilot study, we have 
shown that a 12 -chemokine GES can accurately identify the presence 
of unique, TL-ELNs in metastatic melanoma, which also appear to be 
associated with better patient outcome. We focused on these mela- 
noma cases in particular because such a patient population are con- 
sidered as a candidate to receive immune-based therapies. Moreover, 
this signature has now been interrogated on genomic arrays of 14,492 



SCIENTIFIC REPORT: | 2 : 765 | DOI: 1 0.1 038/srep00765 



3 




C- 



d. 



CD3 



CD4 



CD20 




CD8 - 



CD86 









FoxP3 






100x 




250|jm 








Figure 4 | Representative ectopic lymph node-like structures in melanoma examined by immunohistochemistry. CD20^ B cells are concentrated as a 
follicle, with CD3^, CD4^, and CD8^ T cells appearing in the parafollicular cortex or marginal zones and with some dispersion into the follicle. CD86^ 
antigen presenting cells appear to be dispersed throughout the ectopic lymph node-like structure whereas Fox P3^ staining cells are few in numbers 
(a.,b.,c.). A melanoma metastasis displaying numerous CD20"^ B cell, CD3^ T cell-containing, ectopic lymph node-like structures is also shown (d.). 
Magnification: x200. 
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Table 1 | Pilot Clinical Findings in Melanoma 

Therapy 



Signature" 


TL-ELNs 


Sample No.° 


Tumor Location 


Time (in months) ^ None 


Surgery 


Chemo^ 


XRT 


Other 


+ 


+ 


1 12 


Pancreas 







+ 




Ipi 


+ 


+ 


115 


Lung 


1 14<='d 


+ 






+ 


+ 


116 


Lung 


72c,d 


+ 




+ 




+ 


+ 


113 


Lung 




+ 


+ 






+ 


+ 


118 


Pelvic soft tissue 


65= + 














5 


Liver 


36 




+ 










1 


Brain 


35 


+ 




+ 








3 


Brain 


30 


+ 




+ 


IFN 






6 


Skinf 


13 


+ 












7 


Node 


8 + 











°See Figure 2c. 

'^Time (in months) since primary melanoma to death or present date of evaluation. 
^p=0.008 compared to signature - group (exact Wilcoxon rank sum test). 
'^AliveasofMay 1, 2012. 

^Chemotherapy drugs included: carbo-taxol (1 patient), temozolomide (2 patients), liposomal cytarabine and sagopilone (1 patient each). 
'Subcutaneous metastasis. 

Abbreviations: Ipi, ipilimumab; IFN, interferon-alfa; XRT, external beam radiation therapy. 



distinct primary and metastatic solid tumors (Figure 1), raising the 
intriguing possibility that our findings, earher in CRC and now in 
melanoma, may extend to other histologies as well. The use of the 12- 
chemokine GES was logistically easier, more cost-effective and far 
less time-consuming to predict the presence of TL-ELNs than histo- 
logic analysis based on our experience with CRC. In the earlier study, 
the presence of TL-ELNs often required analysis of numerous stained 
tumor sections, sometimes requiring multiple cuts through the par- 
affin block for positive identification of lymphoid infiltrates consist- 
ent with these structures. Remarkably, this was not the case for the 
melanoma metastases in the current study, which required only a 
single cut section. 

Tumor infiltrating lymphocytes can have profound anti-tumor 
activity in melanoma^'^, so we are in the process of determining the 
functional nature of the TL-ELNs in the laboratory. These structures 
have the general appearance of "typical" peripheral lymph nodes, 
albeit much smaller in size, and are constructed of the necessary 
immune components, including lymphatics (not shown), CD86^ 
antigen presenting cells, and B cell follicles with adjacent T cell 
(CD4^ and CD8^) zones. Of interest, FoxP3^ cells (putative T reg- 
ulatory cells) were found to be very few in numbers. In this regard, 
notably absent from the 12-chemokine GES are CCLl, CCL20, and 
CCL22, which have been shown to selectively recruit and/or main- 
tain T regulatory cells^^"^^. The functional immune nature of the TL- 
ELNs with respect to local (and perhaps also systemic) tumor res- 
ponse and antigen(s) reactivity remains unknown. We plan to 
employ laser capture microdissection with DNA and RNA isolation 
and sequencing methodologies to address questions of clonality and 
functional relationships of the resident B cells and T cells within the 
tumor microenvironment^^'^°, including in depth T cell receptor and 
immunoglobulin transcript repertoire analyses. The presence of 
these unique lymphoid structures also appears to be associated with 
better survival, but this observation will require corroboration in 
larger, prospectively- defined sets of patients. 

In spite of considerable progress made to date, it is well known that 
the majority of melanoma patients treated with immune-based ther- 
apies do not achieve a response or show clinical benefit. We believe 
that our findings raise the intriguing possibility of eventually select- 
ing patients for immunotherapy interventions based on immune- 
related GES from their tumor^\ We plan to evaluate the capacity 
of the 12-chemokine GES to predict the likelihood of patients achiev- 
ing a response of prolonged duration to immune checkpoint inhib- 
itory antibodies (e.g., anti-PD-1 and anti-CTLA-4) and/or cytokines 
(e.g., high-dose IL-2) - first, retrospectively and, later if appropriate. 



prospectively. We will also attempt to extend our findings to other 
solid tumor types (e.g., colon, lung, and ovarian carcinomas), with 
the intent of possibly broadening the use of immunotherapies 
beyond melanoma. Lastly, the use of novel immune- related GES 
may reveal basic information about the mechanism of the anti-tumor 
effects underlying successful immune-based therapies, leading to 
improvements in current immune checkpoint antibody therapy, 
and in the identification and ex vivo expansion of tumor infiltrating 
lymphocytes for broader adoptive immunotherapeutic application. 

Methods 

mRNA microarray analysis. Tumors from patients treated at the Moffitt Cancer 
Center under a University of South Florida Institutional Review Board- approved 
protocol were arrayed on Affymetrix HuRSTA-2a520709 GeneChips (Affymetrix, 
Santa Clara, CA). The chips contain —60,000 probe sets representing —25,037 unique 
genes (Affymetrix HuRSTA-2a520709, GEO: http://www.ncbi.nlm.nih.gov/geo/ 
query/acc.cgi?acc=GPL10379). The expression data for the 14,492 distinct solid 
tumors were normalized using MASS and expression data for the 12 chemokine genes 
(13 probe sets) were extracted. The first principal component, PCI, (explaining 
—53%) was calculated using the Evince™ software (www.umbio.com). The tumor 
samples were grouped by tissue type and the distribution of the PCI score was 
calculated. The expression data for the 120 individual patients, stage IV non- 
locoregional metastasis samples were normalized using RMA and expression data for 
the 12 chemokine genes (13 probe sets) were extracted. A principal component 
analysis (PCA) was performed from these 13 probe sets. The PCI, representing the 
most variability (—57%) within the metastatic melanoma samples, was used to 
represent the chemokine signature. Based on PCI, samples with values >+4 (9 
patient samples) or <— 4 (12 patient samples) were identified as the high and low 
expressers of the chemokine signature, respectively. A heat map was generated using 
Cluster 3.0 with unsupervised, hierarchical method and then visualized using Java 
Treeview on expression data of 12 chemokine genes (13 probe sets) in 21 selected 
samples. The two-sided Student's t test and Bonferroni correction were used to test for 
differences in high (n=9) and low (n= 12) gene expresseion across the genes in the 
chemokine signature. PCA was performed and visualized using Evince™. 

Selection of human tissues. The histological sections corresponding to 21 cases 
(prepared from the mirror-image of the portion of tumor submitted for the mRNA 
microarray analysis) were retrieved from the Moffitt Cancer Center Anatomic 
Pathology Division's repository. All of the specimens were preserved in 10% buffered 
formalin before assessing the presence of microscopically evident host immune 
response. The final pathology report for each case was also reviewed, and the 
pathological data were collected. Linked, annotated but deidentified clinical data (i.e. 
survival) and treatments received were also available in the database. A single, 
random representative section through all selected tissue blocks was evaluated for the 
presence or absence of dense peritumoral infiltrates consistent with TL-ELNs. The 
tissue sections of melanoma metastases were stained using the avidin-biotin complex 
method with retrieval under high pH. Prediluted monoclonal antibodies (mAb) to 
CD3, CD4, CDS, CD20, CD86 and FoxP3 (Ventana Medical Systems, Tucson, AZ) 
were used for the analysis of lymphoid infiltrates. The slides were scanned using the 
Aperio™ ScanScope XT (Vista, CA) with a 200x/0.8NA objective lens at a rate of 2 
minutes per slide via Basler tri-linear-array. These digital slides are housed in the 
Aperio server located in the Moffitt Network Operations Center and are available via 
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the password protected web-based Spectrum database package. This ensures rapid 
quality control of images and analysis by all investigators. Clinical information is 
accessible only to the PI and authorized collaborators and all samples are 
anonymously coded prior to analysis. 

Statistical analysis. Significance of linkages of the 12-chemokine GES profile to 
patient survival was analyzed statistically by the exact Wilcoxon rank sum test. For 
data presented in Figure 1, significance was analyzed by Fisher's exact test. 
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